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ABSTRACT 

We report the results of simultaneous observations of the Vela pulsar in X-rays and radio from the 
RXTE satellite and the Mount Pleasant Radio Observatory in Tasmania. We sought correlations 
between the Vela's X-ray emission and radio arrival times on a pulse by pulse basis. At a confidence 
level of 99.8% we have found significantly higher flux density in Vela's main X-ray peak during radio 
pulses that arrived early. This excess flux shifts to the 'trough' following the 2nd X-ray peak during 
radio pulses that arrive later. Our results suggest that the mechanism producing the radio pulses is 
intimately connected to the mechanism producing X-rays. Current models using resonant absorption 
of radio emission in the outer magnetosphere as a cause of the X-ray emission are explored as a 
possible explanation for the correlation. 

Subject headings: pulsars: X-ray - pulsars: radio - pulsars: individual (Vela) 
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1. INTRODUCTION 

The Vela pulsar (PSR B0833-45) has been well- 
studied. Much observational work has been done 
to understand Vela's emission in individual wave- 
length regions; e.g. the work of Krishnamohan & 
Downs (1983, hereafter KD83) in the radio regime, ob - 
servations by lOegelman. Finlev. fc Zimmermaiil (I1993D 
in the X-ray regime, and studies bv iKanbach et al.l 
lfl99l in the gamma ray regime. Observations of 
Vela's spectrum allow for the p ossibi lity of both po- 
lar cap {D aughcrtv & Harding |1996f ) and outer-gap 
(jCheng. Ruderman. fc Z hang 2000) models of emission. 
Vela's pulse profiles in individual regions have been 
phase-aligned using the phase of the radio pulse across 
the optical, X-ray, and gamma ray wavelength bands 
(Harding et al. 2002, hereafter H02). This article works 
to further relate Vela's high energy emission to its low- 
energy (radio) emission. 

X-ray observations of Vela are challenging. Though 
Vela is the strongest gamma ray source in the sky, 
the pulsar's spectral power drops off in the hard X-ray 
band, making its X-ray emission very difficult to de- 
tect. Additionally, the pulsed emission is overwhelmed 
by the bright but unpulsed background of the X- 
ray emission nebula in which the pulsar is embedded 
(jHelfand. Cotthelf. fc Hal^^[200lh . 

The single-peaked pulse profile of the Vela pulsar in 
radio wavelengths is much simpler than its high en- 
ergy counterparts, although several studies have revealed 
compound emission. KD83 detected peak-intensity de- 
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pendent changes in the pulse-shape with the strongest 
pulses arriving earlier than the averaged profile. They 
conclude that the radio peak is composed of four dif- 
ferent components originating at different heights in the 
emission cone with components lower in the magneto- 
sphere arriving later. In this article one aspect we ex- 
plore is whether a similar connection persists in the X- 
ray regime, i.e. whether the X-ray pulses associated with 
early-arriving radio pulses have a different shape and/or 
a different flux than others. 

Related work on other pulsars includes experiments 
probing the relationship between the Crab pulsar's "gi- 
ant" pulsed emission and its gamma ray emission (Lund- 
gren et al., 1995, hereafter Lu95) or its optical emission 
(Shearer et al., 2003, hereafter Sh03). They reached op- 
posite conclusions. Lu95 observe no correlation within 
their sensitivity, indicating that variations in radio flux 
are caused by changes in radio coherence, which only af- 
fects the radio intensity. Sh03 observe a significant cor- 
relation, and they thus conclude that the increased emis- 
sion in the optical and radio is caused by an increased 
pair production efficiency. 

Patt et al. (1999) study pulse-to-pulse variability in 
the X-ray regime for the Crab pulsar, and they find the 
pulses to be steady to 7%. Using this result, as well 
as previous work showing that the Crab exhibits giant 
radio pulses roughly every two minutes, they conclude 
that the radio and X-ray emission mechanisms are not 
closely related, even though it is likely that the optical 
and X-ray emission regions exist in the same section of 
the mag netosphere (|Patt et al.lll999f ). 

Additional experiments linking pulsar emission in 
the radio and X-ray regimes have been performed 
by Cusumano et al. (2003) and Vivekanand (2001). 
Cusumano et al. show that in PSR B1937+21 there is 
close phase alignment between X-ray pulses and giant 
radio pulses, suggesting a correlation in their emission 
regions. Vivekanand, on the other hand, finds that the 
X-ray flux variations are so much smaller than those at 
radio wavelengths that they are inconsistent with the ex- 
istence of any relationship between the charged emitters 
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in the two wavelength regimes. 

Giant radio pulses have not been shown to exist in 
Vela, but Johnston et al. (2001, hereafter JOl) discovered 
microstructure and 'giant micropulses' in the Vela pulsar. 
The giant micropulses have flux densities no more than 
ten times the mean flux density and have a typical pulse 
width of ^ 400^s. 

By doing a pulse- by-pulse analysis of the Vela pulsar in 
X-ray and radio wavelengths, we will show in this paper 
that the Vela pulsar's X-ray and radio emission mecha- 
nisms must be related. We will discuss the X-ray and 
radio observations in |21 our analysis in the effects of 
scintillation in 521 £^ discussion of interpretations in [JSJ 
and finally our conclusions and related future work in [JB) 

2. OBSERVATIONS 

Our data consist of 74 hours of simultaneous radio and 
X-ray observations taken over three months at the Mount 
Pleasant Radio Observatory in Tasmania and with the 
RXTE satellite. The radio data were acquired during 12 
separate observations between 30 April and 23 August, 
1998. 

The radio data were collected as part of the long term 
monitoring progra m of ten young pulsars, including Vela 
(jLewis et al.l 120 03^. These data were collected with the 
26m antenna at 990.025 MHz using the incoherently de- 
dispersed single pulse system (full description in Dodson, 
McCuUoch, & Lewis, 2002). All individual pulses from 
Vela are detectable, and the pulse height, integrated area, 
and central time of arrival (for the solar-system barycen- 
ter) were calculated from cross-correlation with a high 
signal to noise template in the usual fashion. 

The X-ray data were taken during the same three 
months, yielding 265 ks of usable simultaneous obser- 
vation. For the purposes of this project, only top-layer 
data from RXTE's Proportional Counter Units (PCUs) 
in Good Xenon mode in the energy range of 2-16 keV 
were used. Other filtering parameters included were 
standard RXTE criteria: elevation was greater than 10 
degrees, offset was less than 0.02 degrees, the data were 
taken with at least 3 PCUs on, time since SAA was at 
least 30 minutes, and electronO was less than 0.105. 

3. ANALYSIS 

We filtered the X-ray photon arrival times and trans- 
formed them to the Solar System Ba,rycent er (SSB) using 
the standard FTOOLS (|Blackburnlll995h package. We 
calculated the pulsar phase at the time of each X-ray pho- 
ton, using the radio pulsar-timing program TEMPO^*^, 
and the ephemeris downloadable from Princeton Univer- 
sity We matched each X-ray photon with the radio 
pulse that arrived at the SSB at the same time. The 
precise time span associated with each radio pulse was 
given by our best model for arrival time of the peak of 
the radio pulse, ±0.5 x the instantaneous pulse period 
calculated via the model. Photons arriving on the bor- 
derline were associated with the earlier pulse. We then 
compared pulse profiles for X-rays segregated according 
to the arrival time of the radio pulse. 

Single radio pulses arrive at a range of times around the 
predicted arrival time, as KD83 found. The histogram 
of residual arrival times for radio pulses, relative to the 

See |http: / /pulsar. priiiceton.edu/tempo| 
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Fig. 1. — The number of radio pulses vs phase relative to a 
predictive long-term timing model over the entire data set. The 
dotted lines show the average position of boundaries of the 10 bins 
that were used to make the 10 profiles shown in Figure |2] 

prediction of our best model, is shown in Figure [1] In 
the figure, and in our analysis, the average residual from 
each 5-minute segment of data was subtracted from all 
the data in that segment in order to account for any sys- 
tematic wandering of the pulse arrival times relative to 
our model. The distribution of arrival times is slightly 
skewed, with a tail at late arrival times, so that the mean 
of the distribution is slightly later than the mode (at the 
peak of the histogram), as Figure[l]shows. We divided all 
of the pulses into 10 deciles, by the residual phase of the 
radio pulse, with equal numbers of pulses in each decile. 
Figure [T] shows our division of the residual phase of the 
radio pulse into the deciles. The deciles are well mixed in 
time, i.e. no particular observing epoch dominates any 
decile. Removing the 5-minute average residual elimi- 
nates effects of long-term trends that might appear in- 
dependently in the radio pulse-timing and X-ray photon 
counting data. 

We formed an X-ray pulse profile, integrated from 
2—16 keV, for each of the deciles of radio-pulse arrival 
times, from the X-ray photons associated with each. Fig- 
ure [5] shows the 10 resulting X-ray profiles. Each profile 
contains 13 bins in phase with the radio peak being at 
the left edge of each plot on the border between bins 
13 and 1. The X-ray profiles are significantly different; 
the X-ray pulse changes in shape with the arrival of the 
radio pulse. We denote the ten X-ray profiles by their 
"lateness" : profile 1 comprises the decile of X-ray pho- 
tons from the earliest radio-pulse arrivals, and profile 10 
the decile from the latest radio-pulse arrivals. In partic- 
ular it appears that the first X-ray pulse is sharper and 
stronger in the earlier deciles. Overall 2 distinct X-ray 
peaks are visible in each of the first 5 profiles, whereas 
in profiles 6-10 the two peaks are difficult to distinguish. 
In order to quantify these changes, we performed a num- 
ber of statistical tests on these 10 profiles including a 
full Monte Carlo simulation described at the end of this 
section. 

We observe no significant trend in total X-ray fiux with 
lateness. Figure [3] shows total counts for each of the 10 
deciles. From these data, we determine Pearson's corre- 
lation coefficient of r — —.26 of total X-ray counts with 
radio-pulse lateness with a confidence interval of 54%. 
Nominally, this is not distinguishable from the hypothe- 
sis of zero correlation. Note, however, that Pearson's r 
is not necessarily the best statistic for this comparison, 
as we discuss below. 
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Fig. 2. — Full-period X-ray lightcurves for photons detected during radio pulse arrival times falling in the 10 decile bins shown in Figure 



We do find that the total counts reported in Figure [3] 
are inconsistent with the Gaussian distribution expected 
for this Hmit of the Poisson distribution produced by 
shot noise, at 99.7% confidence, as determined from a 
test. This suggests that there is indeed an evolution 
of X-ray profile with radio lateness. If we consider the 
counts in individual bins, we find that they do deviate 
from the mean more than one would expect from Poission 
statistics, as well. A simple test determines that in 
bins 2 and 3 we can reject the hypothesis of Poisson noise 
around the mean value with 81% and 83% confidence 
respectively. In other words, bins 2 and 3 show larger 
changes than we would expect a priori. The changes 
among the 10 profiles in each of the other 11 lateness 
bins are consistent with Poisson noise. 

The test, however, is not suited to detecting trends. 
In order to detect trends we look to Pearson's correlation 
coefficient, r and its associated confidence interval. Pear- 
son's r has limited validity in our case because it usually 
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Fig. 3. — Total X-ray counts vs 'lateness' decile. 



assumes that both variables are drawn from random dis- 
tributions with nearly Gaussian statistics. In our case, 
lateness is deterministic rather than random and Gaus- 
sian. Shot noise in the profile is random and Gaussian, 
but differences of the X-ray profile among deciles need 
not be. Nonetheless, this is a useful first step to take to 
estimate the significance of any trends. Note also that 
Pearson's r itself, though widely used to measure the 
strength of a correlation, does not judge the existence of 
a correlation. We therefore place more stake in the confi- 
dence interval, although this involv es additional assu mp- 
tions about the variables correlated (jPress et al.lll992[ ) . If 
our visual impression of the profiles is accurate we would 
expect that in the trough, made up of bins 12, 13 and 
1, the correlations are positive (increasing counts with 
increasing lateness). All 3 bins do indeed show positive 
correlations, 0.29, 0.66 and 0.35 respectively with associ- 
ated confidence intervals of 58%, 96% and 68%. In other 
words, 2 of the 3 bins show a significant (greater than 
la — 65%) correlation. Likewise we would expect that 
the first X-ray peak, made up of bins 2 and 3, would show 
negative correlations (decreasing counts with increasing 
lateness). Bins 2 and 3 do in fact show negative correla- 
tions of -0.73 and -0.48 with 98% and 85% confidence, 
consistent with our interpretation of the first peak de- 
creasing with lateness. The only other bin that displays 
a significant correlation is bin 8, which shows a negative 
correlation of —0.63 with 95% confidence. 

Next we looked to see if pairs of bins could be com- 
bined together to better detect the signal. Based on our 
results described above we were curious about whether 
the quantity Yi^a = (ci2 + Cis) - (cii2 + Qia) or Yi^b = 
{ci2 + Qs) — (cii3 + Cii) would show significant trends 
with lateness. Here Q i, Ci2, Qa, q 12, q 13 are the mea- 
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Fig. 4.— 
lateness of the radio pulse. 

sured counts in bins 1, 2, 3, 12, and 13, at lateness i. 
Yi,a and effectively measure, for the 10 profiles, the 
height of the peak minus the trough and are shown in 
Figure U) Via a simple test we can reject with 97% 
confidence the hypothesis that the parent distribution of 
either of the Y^s is a constant at the mean value. More 
importantly Figure [4] shows a clear systematic trend of 
Y with lateness. A line fitted to the data and its as- 
sociated Poisson uncertainty shown in Figure |4] yields a 
slope of m = —297 and —308 respectively with a formal 
1-standard deviation uncertainty of = 76 for both, so 
both represent (not independently) ^ 4 standard devi- 
ation detections of a trend in these data. For the sake 
of completeness we tried all possible differences of pairs 
of adjacent bins. The next highest magnitude slope was 
208, but this and all other significant correlations in- 
cluded some subset of bins 12, 13, 1, 2, and 3. The 
correlation coefficients and associated confidence inter- 
vals for Y, a and Y, b are -0.91 at 99.97% and -0.84 at 
99.74%. 

Given the limited validity of the Pearson's r coeffi- 
cient with non-random, non-normally distributed data, 
we sought a more rigorous estimation of the possibility 
that such a significant trend would arise amongst 10 such 
profiles merely by chance, i.e. merely by random fluctua- 
tion of each bin around its mean. To answer this question 
we performed a Monte Carlo simulation with 10^ realiza- 
tions of the data. Each simulated data set consisted of 
10 profiles, each of 13 bins. The counts in each bin were 
chosen as a Gaussian deviate of the X-ray pulse profile, 
averaged over lateness, and with variance set by Poisson 
statistics. Thus, in the simulated data sets all differences 
among the 10 profiles arose entirely from counting statis- 
tics. In each simulated data set, for every possible pair 
of summed adjacent bins ({1,2 & 3,4}, {1,2 & 4,5} ...{j,k 
& n,o}, with bins {j,k} adjacent and bins {n,o} adjacent: 
66 pairs in all) we computed Yi: 

where i — 1 corresponds to the earliest profile in terms 
of radio phase, and Cij is the number of counts in the jth 
bin of the ith profile. For the set of points {i, Yi} where 
m and a are the slope and uncertainty of the best-fit 
straight line we computed the following statistic: 
r = {m/af 

As for the fit to our data described above, we weighted 
the data by the uncertainties as given by Poisson statis- 
tics. We then compared the largest F for each simulated 



data set to the largest F for the actual data (16.4) and 
found a probability of 0.0024 that the correlation we ob- 
serve could have been obtained by chance. 

We tried other statistics to see if we could find one 
more sensitive to the presence of a correlation like that 
we observed. The most sensitive one included Pearson's 
r as follows: 

F' = + {m/a'f 

where cr' is the ■y/mean(cy + Cik, cu + Cim). F' yielded 
a significance of 0.0007. F' is more sensitive than F to 
the proximity of our data to the fitted line, so it yields 
a smaller probability of false detection. Regardless, we 
retain the more straightforward F as a conservative esti- 
mate of the significance of our result. 

Of use in interpreting our results is knowledge of the 
character of the radio residuals by which we are binning 
the X-ray photons. KD83 did much work in this area, 
but one question they do not address directly is the fol- 
lowing: if a particular pulse is "early" , what is the chance 
that the next pulse will also be early? We calculated the 
autocorrelation function of radio residual, shown in fig- 
ure [51 The function is normalized by the autocorrelation 
at zero lag. The function at a lag of 1 pulse is repre- 
sented by the left edge of the plotted curve, at a value 
of 0.066. The figure shows that the pulsar has very little 
"memory" of the lateness of the previous pulse, although 
it is interesting to note there is finite correlation out to 
40 pulse periods. 

4. SCINTILLATION 

In contrast to effects intrinsic to the pulsar, scintilla- 
tion is unlikely to produce the observed association, be- 
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Fig. 5. — Autocorrelation of radio residual vs lags (pulses). 
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cause it does not affect X-rays; scintillation might erase 
such a correlation but it cannot introduce it. Diffractive 
scintillation for Vela at our frequency has a characteristic 
bandwidth of 2 kHz (Gupta 1995), far smaller than our 
observing bandwidth of 6.4 MHz. We therefore average 
over ~3200 independent scintillation elements. Refrac- 
tive scintillation modulates flux densit y over a wide band - 
width and has timescale ~25 days ( Desai et al.lll992l ). 
This was shorter than the total span of the observations, 
but much longer than the span of any one observation. 
To combat any effect that this might have on the ob- 
served arrival times, we defined the 10 deciles separately 
for each observation; i.e. the specific values of residual 
that separated each of the 10 bins were calculated for 
each individual radio observation. Again, these cutoffs 
were defined in such a way that an equal number of ra- 
dio pulses was associated with each decile. 

We assumed that the pulsar was stationary relative 
to a refractive scintillation element during each single 
observation, and further tests to ensure that the length 
of the observation was not a factor were performed by 
renormalizing the cutoffs using both one hour and five 
minute timespans. Reanalyzing the X-ray data in one 
hour segments did not significantly change our results, 
and the five minute spans were found to be too short to 
accurately represent Vela's emission. 

5. DISCUSSION 

The early radio emission could result from coherent 
radiation along a different set of field lines (i.e. more 
leading) or from radiation at a higher altitude, or both. 
The work of KD83 suggests both. The early radio emis- 
sion may also result from stochastic fluctuations in the 
radio beam intensity which would lead to pulse arrival 
time changes by changing the shape of the radio beam 
as it takes a finite time to sweep across our line of sight. 

Petrova (2003, and references therein) offers a physical 
model that could explain the observed relationship be- 
tween the radio and X-ray emission. Her model suggests 
that resonant absorption of radio emission from the outer 
magnetosphere leads to an increase in the pitch angles 
and momenta of the secondary pairs, which then leads 
to optical and higher energy emission by spontaneous 
synchrotron radiation. How could this model produce a 
changing X-ray pulse shape? Due to the effects of ro- 
tation of the magnetosphere and aberration, the early 
radio emission can cross a larger number of higher alti- 
tude field lines and at larger angles, thereby maximizing 
the opportunity for absorption by particles on those field 
lines, and therefore the production of high-energy radi- 
ation. Conversely a late radio photon on a path almost 
directly along the magnetic pole may escape the magne- 
tosphere with many fewer interactions, since it will cross 
fewer open field lines and at smaller angle. The details of 
the above are somewhat unimportant as our knowledge 
of the magnetosphere and the plasma therein is limited, 
but the point is that as different parts of the radio beam 
are active, resonant absorption may happen at different 
rates in different parts of the magnetosphere, causing 
continuous change in the shape of the observed X-ray 



emission. 

More generally our results imply a connection between 
the radio and X-ray emission mechanisms for Vela that 
is not consistent with outer gap models. In these models, 
the high energy emission results from a gap connection to 
the pole opposite from that producing the radio emission. 
The pole and outer gap associated with the same set of 
field lines are not visible to one observer. It is not clear 
how a correlation could exist between the radio and high 
energy regimes in these models. 

The giant micropulse emission observed by JOl would 
cause a single radio pulse to arrive about 1 ms early, 
so it is realistic to consider the possibility that the gi- 
ant micropulse emission is primarily responsible for the 
early arrival of the radio pulse. However, out of 20,085 
pulses that JOl observed at 1413 MHz, 14 of them con- 
tained giant micropulses. So the giant micropulses may 
be influencing the flrst of our 10 deciles, but cannot con- 
tribute to the effect observed in the other 9 deciles. We 
conclude that Vela's giant micropulse emission cannot be 
responsible for the effect presented here. 

6. CONCLUSIONS AND FUTURE WORK 

We have evidence that links features of Vela's X-ray 
emission with features of its radio emission. First, we 
find that X-ray pulses associated with early radio pulses 
show stronger emission at the main X-ray peak which 
is the sharper of the two. Similarly X-ray pulses associ- 
ated with later radio pulses show stronger emission at the 
trough following the 2nd X-ray peak, a region in phase 
near the radio peak. The trend we measure has a 0.2% 
probability of appearing in the data by chance. We con- 
clude that there is a close relationship between X-ray and 
radio emission in the Vela pulsar. 

We plan to further characterize the relationship be- 
tween the radio and high energy emissions of pulsars 
to identify their origins and constrain magnetospheric 
models. In particular, we will explore the dependence of 
radio-to-X-ray correlations on the radio frequency and 
polarization properties of individual Vela pulses, both of 
which carry information about emission altitudes. Sim- 
ilar observations of other pulsars also promise useful in- 
sights as probes of different magnetic field strengths and 
emission/ viewing geometries. 
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